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J^* Abstract 

This fluid dynamics video shows the breakup of a droplet in a 
stationary homogeneous and isotropic turbulent flow. We consider 
t-H droplets with the same density of the transporting fluid. The droplets 

^ and the fluid are numerically modelled by means of a multicompo- 

nent Lattice-Boltzmann method. The turbulent fluid is maintained 
through a large scale stirring force and the radius of stable droplets, 
for the parameters in our simulation, is larger than the Kolmogorov 
scale. Events of droplet deformation, break-up and aggregation are 
clearly visible from the movie. With the present database droplet evo- 
lution can be studied from both an Eulerian and Lagrangian point 
t—H of view. The Kolmogorov-Hinze criteria for droplets break-up can be 

tested also by means of simulations with different viscosity contrast 
• t-h between the two components. 

1 Introduction 
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Droplet emulsions are key to many natural and industrial processes. In pres- 
ence of an external flow, droplets undergo deformation, breakup, and coagu- 
lation. In a turbulent flow, breakup of droplets larger than the Kolmogorov 
scale is governed by the interplay between surface tension and turbulent 
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Figure 1: The initially spherical droplet (left panel) get strechted (central 
panel) by switching on the stirring force and generation of a turbulent flow. 
The righ panel shows the initial droplet immediately after its first breakup. 

pressure fluctuations. Balancing surface tension and pressure fluctuations 
allows to estimate the minimum unstable droplet diameter that can under- 
goes breakup [U [2] . Turbulent pressure-driven and surface-tension stresses 
across a droplet of typical diameter d and with surface tension a can be 
estimated as: 

Tturb oc pc^Sdu) 2 ); r tens oc a/d 

where with p c we indicate the density of the continous phase and with 
((Sdu) 2 ) we indicate a typical velocity fluctuations across the droplet. The 
ratio between the two gives the Weber number estimated on the droplet 
diameter: 

We(d) = dpc(( ^ )2) . 

a 

The maximum stable mean droplet diameter can be estimated by assuming 
K41 turbulent statistics for velocity increments: ((8du) 2 ) ~ £ 2 / 3 d 2 / 3 where e 
is the energy dissipation. From the condition of Weber number order unity, 
one obtains: 

/ rr \ 3 / 5 
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Corrections to finite viscosity of the dispersed phase can also be added [3j. 
The present simulation allows to study droplets breakup conditions e.g. test- 
ing the Kolmogorov-Hinze criteria at changing the Reynolds number and 
the viscosity ratio between the dispersed and continuos phase. Specifically 
it is also possible to study the fluctuations around the mean average stable 
droplet radius ([!]) induced by intermittent fluctuations of the velocity field. 

In the movie various stages of droplet breakup process are shown. The 
turbulent flow is homogeneous and isotropic and it is numerically integrated 
by means of a multi-component Lattice-Boltzmann (MCLB) simulation in- 
tegrated in order to evolve consistently droplets and the advecting fluid [I]. 
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Figure 2: Different phases in the time evolution of a large droplet in a 
turbulent flow (•). Initially the droplet deforms until it is broken it two 
droplets. This event is rapidly followed by several other multiple fragmen- 
tation events until a stationary droplet dispersion is attained. In the final 
stages the droplets continuously undergo deformation, break-up, and coag- 
ulation. Here the time is non dimensionalized by the large eddy turnover 
time r eddy . 



The multicomponent flow is simulated using a standard Shan-Chen model 
[5], the turbulence is generated by applying a large scale forcing, and peri- 
odic boundary conditions are used. The parameter used in our simulation 
are given in table [l} 

We start the simulation with the fluid at rest and all the droplet volume 
fraction residing into a large, single spherical droplet. From the beginning 
of simulation, the large scale stirring start to be applied and the turbulence 
level grows continuosly during roughly one large scale eddy turnover time. 
In the initial phase the large drop gets mildly deformed, then stretched, 
and finally broken into smaller droplets (see Figure [I]). At this point the 
fragmentations process continues for about another 1 2 eddy turn over 
times. After this period the distribution of droplets is almost stationary 
and is characterized by continuos coagulation, deformation and breakups 
events due to the continuous presence of underlying turbulent fluctuations. 
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N a v 




d 


d c 


d lbm 


512 1.6 x 10" 3 5 x 10" 3 


8.5 x 10" 3 29 


50 


24.16 


25 ±2 



Table 1: Parameters of the numerical simulation from which the movie was 
produced. N is the number of grid points along each direction, a is the 
surface tension, R\ = u rms \/v is the Taylor-microscale Reynolds number, 
A = y/Eyn is the Taylor-microscale, E is the average fluid kinetic energy 
and £1 is the average enstrophy. do is the initial droplet diameter, d c is the 
critical droplet diameter obtained from the theoretical estimate (see text) 
and d^ BM is the critical droplet diameter as measured from our simulation. 
The Reynolds number was kept low on purpose in order to keep the dis- 
sipative scale relatively large: this allow us to have an interface thickness 
smaller than the viscous scale. 

For very long simulation one may experience some droplets evaporation. The 
different phases in the time evolution of droplets are schematically illustrated 
in Fig. [2} 

Interesting questions concern the rate of collision and the nature of defor- 
mation and breakup events. It is also of interest to quantify the evolution 
of the statistical properties of droplets diameter and the correlations be- 
tween stresses on the droplet and break-up dynamics. These results will be 
reported in a forthcoming article [B] . 

We acknowledge support from the Juelich supercomputing center (Ger- 
many) for the computational resources. 
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